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H2B monoubiquitylation [H2Bubl}, which is required for multiple methylations of both H3K4 and H3K79, has been 
implicated in gene expression in numerous organisms ranging from yeast to human. However, the molecular crosstaIi< 
between H2Bubl and other modifications, especially the methylations of H3K4 and H3K79, remains unclear in verte- 
brates. To better understand the functional role of H2Bubl, we measured genome-wide histone modification patterns in 
human cells. Our results suggest that H2Bubl has dual roles, one that is H3 methylation dependent, and another that is H3 
methylation independent. First, H2Bubl is a 5'-enriched active transcription mark and co-occupies with H3K79 methylations 
in actively transcribed regions. Second, this study shows for the first time that H2Bubl plays a histone H3 methylations- 
independent role in chromatin architecture. Furthermore, the results of this work indicate that H2Bubl is largely positioned 
at the exon-intron boundaries of highly expressed exons, and it demonstrates increased occupancy in skipped exons 
compared with flanking exons in the human and mouse genomes. Our findings collectively suggest that a potentiating 
mechanism links H2Bubl to both H3K79 methylations in actively transcribed regions and the exon-intron structure of 
highly expressed exons via the regulation of nucleosome dynamics during transcription elongation. 

[Supplemental material is available for this article.] 



Post-translational modification of the four core histories is a ubiq- 
uitous process during the regulation of gene activation and re- 
pression. Histone modifications are also involved in various cel- 
lular processes, including alternative splicing (Kolasinska-Zwierz 
et al. 2009; Luco et al. 2010). Several histone modifications can 
engage in crosstalk with other histone modifications; in either cis 
or trans. One well-known example of histone crosstalk in trans 
is the requirement of H2B C terminus monoubiquitylation (H2B 
monoubiquitylation, H2Bubl) for the multiple methylations 
of H3K4 and H3K79 (Dover et al. 2002; Sun and Allis 2002; 
Shahbazian et al. 2005), which is evolutionarily conserved from 
yeast to human. In Saccharomyces cerevisiae, monoubiquitylation 
of lysine residue 123 in histone H2B was shown to play a positive 
role in H3K4 and H3K79 methylations, which are involved in 
transcription initiation and elongation, respectively (Briggs et al. 
2002; Sun and Allis 2002). In mammals, H2Bubl coupled to RAD6- 
mediated transcription has been reported to directly stimulate 
H3K4 methylation (Kim et al. 2009), and synthetically ubiquitylated 
H2B stimulates DOTiL-mediated H3K79 methylations in vitro 
(McGinty et al. 2008; Kim et al. 2009; Oh et al. 2010). Despite many 
efforts to infer the roles of histone modifications, however, the 
crosstalk between H2Bubl and the multiple methylations of H3K4 
and H3K79 is not fully understood, especially in mammals. 
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The signaling between H2Bubl and histone H3 methylation 
is clearly linked to chromatin dynamics during transcription elon- 
gation, but the manner with which H2Bubl -containing nucleo- 
somes regulate chromatin architecture independently of histone 
H3 methylation during transcription elongation is unclear. In this 
regard, Fleming et al. (2008) found that H2Bubl plays specific roles 
in GALl/10 gene expression in budding yeast, and these roles are 
independent of its regulation of H3 methylation. In addition, Pavri 
et al. (2006) demonstrated that H2Bubl directly stimulates FACT 
(FAcilitates Chromatin Transcription) -mediated transcription us- 
ing an in vitro chromatin-based transcription system. The ubiq- 
uitin moiety of histone H2B is removed by the deubiquitylases, 
UBP8 and UBPl 0 (Weake and Workman 2008). Loss of the H2B E2/ 
E3 ubiquitin ligase RAD6/BRE1 or substitution of the H2B C-ter- 
minal lysine residue with arginine leads to transcription defects, 
suggesting that the dynamic turnover of H2Bubl is critical to both 
transcription and gene silencing. Recent studies have proposed 
that H2Bubl regulates chromatin dynamics by enhancing nucle- 
osome stability (Chandrasekharan et al. 2009) and disrupts local 
and higher-order chromatin compaction (Fierz et al. 2011). These 
findings strongly suggest that H2Bubl is involved in additional 
regulatory pathways beside frans-crosstalk with the methylations 
of histone H3. To understand the genome-wide characteristics of 
H2Bubl in mammals, Minsky et al. (2008) developed a mammalian 
H2Bubl -specific monoclonal antibody. Subsequent use of ChlP- 
chip and ChlP-seq analysis showed that H2Bubl is associated 
with the transcribed regions of highly expressed genes in HeLa 
cells (Shema et al. 2008). However, the prior data are insufficient 
for an investigation of the genome-wide properties and unique roles 
of H2Bubl in cooperation with other modifications. Therefore, we 
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herein generated high-throughput sequencing (ChlP-seq) data for 
the genome- wide occupancies of H2Bubl, nucleosome, H3K4me3, 
H3K36me3, H3K79mel/2/3, H3ac, and mRNA-seq in human em- 
bryonic carcinoma (NCCIT) cells. Furthermore, to investigate the 
effects of H2Bubl reduction, we performed mRNA-seq and ChlP-seq 
of H2Bubl in NCCIT cells transfected with i^NF20-small interfering 
RNA (siRNF20), which decreased the H2Bubl signal. We found that 
H2Bubl is a 5 '-enriched active transcription mark and that the- 
methylations of H3K79 co-occupied with H2Bubl in highly 
expressed genes. Notably, H2Bubl was strongly enriched at the 
exon-intron boundaries of highly expressed exons and appeared 
to correlate with exon skipping. 

Results 

H2Bubl is a 5'-enriched active transcription mark 

To determine genome-wide patterns of histone modifications, in- 
cluding histone H2B monoubiquitylation; we carried out ChlP-seq 
analysis using undifferentiated human teratocarcinoma NCCIT 
cells. NCCIT cells were established from a mediastinal mixed germ 
cell tumor (Teshima et al. 1988) at an intermediate stage between 
seminoma and embryonal carcinoma (Damjanov 1993). NCCIT is 
a developmentally pluripotent cell line that can differentiate into 
derivatives of all three embryonic germ layers and extra-embryonic 
cell lineages (Damjanov et al. 1993). H2Bubl demonstrated rela- 
tively strong signals near the 5' regions of actively transcribed 
genes, and then gradually decreased along the transcribed region; 
this is shown for the UBE4B gene (ubiquitin conjugation factor E4 
B) in Figure lA, as an example. Similar patterns were also observed 
for H3K79 methylations, while H3K4me3 and H3ac were enriched 
in the regions surrounding the transcriptional start sites (TSS). 
These patterns, as well as all other observed modification patterns, 
were consistent with previously reported results in human cells 
(Barski et al. 2007; Wang et al. 2008). 

To infer the genomic regions that showed H2Bubl enrichment, 
we detected H2Bubl peaks using the MACS peak detection software 
(see Supplemental Methods). In agreement with the 5 '-enriched 
pattern we had observed for H2Bubl, most of the H2Bubl peaks 
were located in the 5' UTR and open reading frame (ORF) regions 
(Fig. IB). This result suggests that the specific function of H2Bubl 
is localized to transcribed regions, not promoter regions, further 
indicating that H2Bubl may play a role in potentiating transcrip- 
tional elongation. The H2Bubl peaks highly overlapped with the 
H3K79 methylation peaks; -24% of the H3K79me3 peaks and 
H2Bubl peaks matched, while H3ac and H3K4me3 had similar 
numbers of peaks (Fig. IC). Gene set analysis using DAVID in- 
dicated that genes doubly targeted by both H2Bubl and H3K79 
methylation peaks were enriched in protein modification, me- 
tabolism, and mRNA-related biological functions (Supplemental 
Fig. lA). 

Since previous studies have suggested that H2Bubl is associ- 
ated with highly expressed genes (Minsky et al. 2008; Shema 
et al. 2008), we investigated whether H2Bubl is a transcription- 
dependent active mark in human cells. H2Bubl demonstrated 
strong signals in highly expressed genes (Fig. ID, black line), and 
these signals decreased along the transcribed regions (Fig. IE). The 
H2Bubl signal showed a more than twofold enrichment in highly 
expressed genes compared with lower-expressed genes. To examine 
whether the strong H2Bubl signal was associated with increased 
gene activation, we performed mRNA-seq for siRNF20-transfected 
NCCIT cells. Our results revealed that H2Bubl, H3K4me3, and 



H3K79me3 levels decreased in the RNF20 knockdown cells 
(Supplemental Fig. 2). Since the H3K4me3 and H3K79 methyla- 
tions, which are known to function as gene activation signals, 
decreased with the reduction of H2Bubl signal, we hypothesized 
that the highly expressed genes targeted by H2Bubl were expected 
to demonstrate down-regulation. However, using a twofold ex- 
pressional change as a cutoff value, we found that although 4654 
genes were down-regulated in siRNF20-transfected NCCIT cells, 
a comparable number (5304) were up-regulated (Fig. IF). 

To investigate the biological functions of specifically regu- 
lated genes in siRNF20-transfected NCCIT cells, we conducted 
gene set analysis (Supplemental Fig. IB). According to the GO 
(Gene Ontology), the down-regulated genes were enriched for cell 
cycle-related functions related to tumor suppression. This result is 
consistent with the previous suggestion that H2Bubl may func- 
tion as a tumor suppressor (Shema et al. 2008). The up-regulated 
genes were also enriched for developmental regulation-related genes, 
suggesting a new possible role for H2Bubl in developmental pro- 
cesses. The number of up-regulated genes observed in siRNF20- 
transfected NCCIT cells also raises the possibility that H2Bubl may 
have an additional role in gene regulation beyond its crosstalk with 
H3 lysine methylations. However, although many genes were 
up-regulated in siRNF20-transfected NCCIT cells, the H2Bubl- 
targeted genes tended to be significantly decreased (Fig. IG, 
brown and orange), while those not targeted by H2Bubl tended 
to be slightly increased (Fig. IG, yellow-green and green). The pre- 
sence of gene expression-dependent H2Bubl patterns and more 
significant changes in H2Bubl -targeted genes in the siRNF20- 
transfected NCCIT cells suggest that H2Bubl is a 5 '-enriched gene 
activation mark. 

rra/75-crosstalk between H2Bubl and H3K79 methylations 
is associated with gene expression 

H2Bubl has demonstrated frans-crosstalk with the methylations of 
H3K4 and K79 in a variety of organisms ranging from yeast to hu- 
man (Dover et al. 2002; Sun and Allis 2002; Kim et al. 2005, 2009; 
Shahbazian et al. 2005). We verified these relationships in NCCIT 
cells (Supplemental Fig. 2), and then set out to investigate whether 
this trans-crosstalk generated genome- wide co-occupation patterns 
and/or was associated with gene activation. We first clustered the 
H2Bubl patterns in the regions surrounding the TSS (including 
5 kb upstream/downstream from the TSS) using /c-means clustering 
(see Supplemental Methods). We found three distinct patterns 
(Fig. 2A). After sorting genes according to the H2Bubl -clustered 
results in NCCIT cells, we compared our results with those for 
H2Bubl in HeLa cells (Shema et al. 2008) as well as those for 
H3K79mel/2/3, H3K4me3, and H3K36me3 in NCCIT cells (Fig. 
2A). The H2Bubl clustering patterns were fairly consistent between 
NCCIT and HeLa cells, indicating that these cluster patterns are 
unlikely to be cell-type specific. The methylations of H3K79 dem- 
onstrated genome-wide correlation patterns with H2Bubl. How- 
ever, although H2Bubl is a prerequisite for H3K4me3 (Kim et al. 
2009), H3K4me3 did not highly overlap with H2Bubl. The tri- 
methylation patterns of H3K36 (used here as a control) were in- 
dependent of H2Bubl. These results were consistent with the low 
number of overlapping peaks between H2Bubl and H3K4me3/ 
H3K36me3 (Fig. IC). The discrepancy between H3K4me3 and 
H3K79mel/2/3 suggests that H3K4me3 is affected by other path- 
ways in addition to an H2Bubl -dependent pathway; these might 
include crosstalk with the lysine acetylations of H3 (Hung et al. 
2009), counter-correlation with H3R2 methylation (Guccione et al. 
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Figure 1. Varying genonnic properties of the H2Bub1 -enriched regions. (/4)The ChlP-seq results for t76f4S(ubiquitin conjugation factor E4B) are shown. 
(K-axis) Read counts of the sequence tags in 30-bp increments. H2Bub1 was enriched at the 5' end and gradually decreased along the transcribed region, 
similar to the methylation pattern of H3K79. (B) MACS analysis recognized 14,679 H2Bubl peaks and assigned them to 41 73 genes. (C) The H2Bubl 
peaks overlapped with the H3K79me3 peaks (24% of the H3K79me3 peaks) in 2765 genes. (D,E) H2Bubl signals in the regions surrounding the tran- 
scription start site (TSS) and transcription end site (TES). (/-axis) H2Bubl signal; (x-axis) relative position (bp). Nucleosomes were well positioned in the 
regions surrounding the TSS, while H2Bubl signals were depleted in these regions. The gene expression levels were divided into three groups: highly 
expressed (black lines, top 30%), lower expressed (green lines, bottom 30%), and intermediately expressed (red lines, middle 40%). (f) Comparison of 
gene expression changes between wild-type and siRNF20-transfected NCCIT cells. (Red and blue) Genes demonstrating twofold or more up-regulation 
and down-regulation, respectively. (C) Gene expression changes among the H2Bubl -targeted genes. Among —25,000 genes, 3824 were H2Bubl- 
targeted genes, and the remaining 21,310 genes were defined as H2Bubl off-targeted genes. P-values were calculated using the Kolmogorov-Smirnov 
(KS) test (Stephens 1970). Parentheses indicate median values. 



2007), and several other related pathways. Clearly, the bidirectional 
pattern of H3K4me3 in the regions surrounding the TSS and the 
broader distribution of H2Bubl in transcribed regions provide 
strong evidence that there is a discrepancy between H2Bubl and 
H3K4me3. However, the hDOTlL-mediated H3K79 methylations 



are tightly regulated by H2Bubl in the same nucleosome (McGinty 
et al. 2008). In agreement with these clustering patterns, the H3K79 
methylations and H2Bubl, but not H3K4me3, showed strongly 
correlated patterns of average modification signals in the gene body 
regions (Supplemental Fig. 3A-E). 
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Figure 2. Genome-wide correlation patterns between H2Bubl and H3K79 methylations. (A) Three distinct H2Bubl patterns were identified for highly 
expressed genes (top 30%, 7540) in the regions surrounding the TSS, as assessed by /c-means clustering. Other modification patterns are presented in 
the gene order shown for H2Bub1 . Cluster A shows strongly enriched bidirectional patterns, while both clusters B and C were enriched downstream from 
the TSS. (B) The correlation patterns between H3K79me2 and H2Bubl affected gene activation. Genes were divided into five groups according to their 
gene expression levels (purple: highly expressed genes; yellow: lower expressed genes). Pearson correlation coefficients were calculated between H3K79 
methylations and H2Bub1 within each gene, (/-axis) Average correlation coefficient of each group. (C) Boxplots of gene expression levels in FPKM scores 
for different combinations of H2Bubl - and H3K79-methylation-targeted genes. Modification-targeted genes were defined as those containing at least one 
corresponding modification peak. (Values in parentheses) Number of genes. (D) Gene set analysis for all combinations of H2Bubl- and H3K79 meth- 
ylation-targeted genes using DAVID. PANTHER biological processes are presented with corrected P-values. 



Since H2Bubl demonstrated genome-wide correlations with 
the methylations of H3K79, which are regarded as gene activation 
marks, we investigated whether this genome-wide correlation 
pattern was associated with changes in gene activation. If trans- 
crosstalk between H2Bubl and H3K79 methylations is important 
in gene activation, H3K79 methylations should show strong cor- 
relations with H2Bubl in highly expressed genes. Thus, gene ex- 
pression levels and correlation coefficients were compared be- 
tween H2Bubl and the methylations of H3K79 in the gene body 
regions. Among the mono-, di-, and trimethylations of H3K79, 
dimethylation showed the most significant relationship, with highly 
expressed genes (purple, top 20%) demonstrating more than twofold 
higher correlation coefficients with these genes compared with the 
lower-expressed genes (Fig. 2B). This tendency was maintained for 
the mono- and trimethylations, suggesting that frans-cros stalk 
between H2Bubl and the methylations of H3K79 is related to 
gene activation (Supplemental Fig. 3F,G). 

To test the effect of frans-crosstalk on gene activation, we ex- 
tended our analysis to genes that are cotargeted by H2Bubl and the 
methylations of H3K79. Using the modification peaks defined in 
Figure 1, B and C, we identified genes that were solely targeted by 
H2Bubl, those that were doubly targeted (H2Bubl-H3K79mel, 
H2Bubl-H3K79me2, and H2Bubl-H3K79me3), and those that 
were simultaneously targeted (H2Bubl-H3K79me 1/2/3). We found 
that the H2Bubl -targeted genes tended to be simultaneously cotar- 
geted by H3K79 mono-, di-, and trimethylations (around 60% of the 
4173 H2Bubl -targeted genes). Furthermore, the genes that were 



simultaneously targeted (H2Bubl-H3K79mel/2/3) were highly 
activated compared with H2Bubl solely or doubly targeted genes 
(Fig. 2C). In agreement with this result, high signals for both 
H2Bubl and H3K79 methylations in the gene body regions were 
associated with high gene expression levels (red spots in Supple- 
mental Fig. 3C-E). 

Despite this apparent frans-crosstalk between H2Bubl and 
H3K79 methylations, however, it was clear that genes were tar- 
geted by different combinations of these modifications. In this 
regard, we wondered what properties might be associated with 
different combinatorial patterns of these modifications. All 16 
possible gene groups were generated according to the presence of 
each modification and gene set analysis was conducted. Each 
combinatorial pattern was characterized by distinct biological 
functions (Fig. 2D). 

H2Bubl is highly enriched at the exon-intron boundaries 
of highly expressed exons 

Recent yeast studies showed that C-terminal mutants of H2B 
(K123R in budding yeast and Kl 19R in fission yeast) demonstrated 
more severe phenotypes and modulated transcription among more 
genes compared with H3K4R or H3K79R mutants, indicating that 
H2Bubl plays one or more roles in transcription elongation inde- 
pendent of its regulation of histone H3 methylation in yeast 
(Tanny et al. 2007; Fleming et al. 2008). Since several papers have 
emphasized the role of H2Bubl in modulating nucleosome ar- 
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chitecture (Chandrasekharan et al. 2009; Batta et al. 2011; Fierz 
et al. 2011; Schulze et al. 2011), we herein sought to identify cor- 
relation patterns between H2Bubl and nucleosome architecture. 

First, we focused on the relationship between H2Bubl signals 
and exon-intron structure. To investigate the various modification 
patterns in the regions surrounding the exon-intron boundaries, 
we used /c-means clustering (see Supplemental Methods). As shown 
in Figure 3A, four distinct patterns were clearly identified in the 
exon boundary regions (encompassing 500 bp upstream of and 
downstream from the exon 5' end): high density in the central 
region (cluster A), slightly biased to the exonic region (cluster B), 
slightly biased to the intronic region (cluster C), and low density in 
the central region (cluster D). 

This property of H2Bubl was also observed in HeLa cells 
(Supplemental Fig. 4A). Because exon boundaries are also nucleo- 
some enriched (Andersson et al. 2009; Tilgner et al. 2009), H2Bubl 
density was normalized with respect to nucleosome occupancy 
(see Supplemental Methods). H2Bubl was still enriched at the 
exon boundaries after this normalization (Fig. 3 A, H2Bubl/Nuc), 
implying that this enrichment was not attributable to nucleosome 
occupancy. Unlike H2Bubl, no patterns were recognized in H3K4me3 
and H3K79 methylations (Fig. 3A; Supplemental Fig. 4A), dem- 
onstrating that the enrichment of H2Bubl at the exon boundaries 
was independent of crosstalk with H3 methylations. We also con- 
firmed that H3K36me3 and nucleosome were enriched at the exon 
boundaries (Fig. 3A) as previously reported (Kolasinska-Zwierz et al. 
2009; Schwartz and Ast 2010). 

Next, we examined whether the degree of H2Bubl enrich- 
ment depended on the expression levels of exons. Indeed, highly 



expressed exons showed stronger enrichment at both the 5' and 3' 
ends of the exon boundaries (Fig. 3B,C, black lines). However, in 
intermediately or lower-expressed exons, H2Bubl demonstrated 
only very slight enrichment or no specific patterns at the exon 
boundaries (Fig. 3B,C, red and green lines). Similar results were 
obtained in HeLa cells (Supplemental Fig. 4B,C) and mouse myo- 
blast cells (Vethantham et al. 2012; Supplemental Fig. 4D,E). In- 
terestingly, Pearson correlation coefficient showed a slight positive 
correlation between the H2Bubl signals and exon expression levels 
(r = 0.05, P = 4.2 X 10"^^). However, this was less of a correlation 
than that of H3K36me3 (r= 0.19, P= 0) (data not shown), suggesting 
that H2Bubl is important in highly expressed exons, but that it is 
less prominent than H3K36me3 in determining exon expression 
levels. 

It is possible that the high H2Bubl densities at the boundaries 
of highly expressed exons are just the result of the overall higher 
level of H2Bubl throughout the bodies of highly expressed genes 
compared with intermediately or lower expressed genes (Fig. 3D). 
However, regardless of the gene expression levels, highly expressed 
exons showed stronger H2Bubl signals (Fig. 3D, High) compared 
with the remaining exons (Fig. 3D, Others). This tendency was 
more significant in the intermediately or lower-expressed genes 
(Fig. 3D, blue bars), supporting the notion that H2Bubl enrich- 
ment plays a unique potentiating role in highly expressed exonic 
regions. 

To validate the role of H2Bubl in highly expressed exons, we 
used wild-type and siRNF20-transfected NCCIT cells, and com- 
pared the expression levels among exons with high H2Bubl sig- 
nals (top 30% of exons according to H2Bubl signal) and those with 
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low H2Bubl signals (bottom 30% of exons). The expression 
levels of exons with high H2Bubl signals significantly decreased 
when H2Bubl expression was reduced by siRNFZO transfection 
(Fig. 3E, brown and orange). In contrast, the expression levels 
of exons with low H2Bubl signals decreased less significantly 
in siRNFZO-trasnfected NCCIT cells (Fig. 3E, yellow-green and 
green). Thus, our findings collectively suggest that H2Bubl is 
strongly enriched at the exon boundaries of highly expressed 
exons and that this property is evolutionarily well conserved in 
mammals. 

We next examined whether distinct H2Bubl patterns at the 
exon boundaries (Fig. 3A) could be associated with exon-intron 
structure. When patterns of exon length, gene strand, upstream 
intron length, and exon position were analyzed according to clus- 
ters, only upstream/downstream intron length and exon position 
showed cluster-specific patterns. Exons located following short 
upstream/downstream introns tended to be included in cluster A 
or D compared with cluster B and C (Supplemental Fig. 4F,G). 
Additionally, more exons near the 3' end tended to be assigned to 
cluster D (Supplemental Fig. 4H). This tendency may be associ- 
ated with the 5 '-bias of the H2Bubl pattern. Most of the exons in 
cluster D were marked by H3K36me3 (84% of the exons in cluster 
D were included in clusters A-C of H3K36me3). These data in- 
dicate that H3K36me3 may play a complementary role in exons 
unmarked by H2Bubl (Supplemental Fig. 41). 

H2Bubl is highly enriched in skipped exons in mammals 

Previous studies reported that two exon determination marks, 
nucleosome and H3K36me3, are less prominent in skipped exons 



(Hon et al. 2009; Choi 2010), and an association between histone 
modifications and alternative splicing has been experimentally 
validated (Luco et al. 2010). Since H2Bubl is also enriched at exon 
boundaries, we compared H2Bubl signals between skipped exons 
and flanking exons (see Methods) to determine whether H2Bubl 
was correlated with skipped exons. Interestingly, H2Bubl signals 
were higher for skipped exons in human cells (Fig. 4A,B). Consis- 
tent with previous results, reduced nucleosome and H3K36me3 
signals were also observed in the skipped exons (Fig. 4C; Supple- 
mental Fig. 5 A), whereas no specific pattern was observed for the 
control, H3K79me2 (Fig. 4D), or H3K4me3/H3K79mel/H3K79me3 
(data not shown). To test whether this result was specific to human 
cells, genome-wide H2Bubl data from mouse and Drosophila cells 
were also analyzed. Notably, an evolutionarily conserved pattern 
was observed in the mouse (Fig. 4E), while an opposing result (i.e., 
reduced signals in skipped exons) was found in Drosophila cells 
(Fig. 4F). Exon skipping was also correlated with the clustered pat- 
terns of H2Bubl (Fig. 3A), as skipped exons were slightly enriched 
in clusters A and B compared with clusters C and D (Supple- 
mental Fig. 5B). Interestingly, the tendency of H2Bubl speci- 
ficity in skipped exons was stronger in clusters A and B compared 
with clusters C and D (Supplemental Fig. 5C). These cluster-de- 
pendent patterns support a putative H2Bubl role in exon skip- 
ping, since clusters A and B include exons enriched by H2Bubl in 
their exonic regions. 

To validate the association between H2Bubl and exon skip- 
ping, the number of skipped exons was measured in siRNF20- 
transfected NCCIT cells. Of the 34,763 exons skipped in wild-type 
cells, 21,804 were persistently skipped in siRNF20-transfected NCCIT 
cells, while the remaining 12,959 exons were included. Thus, the 
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Figure 4. Skipped exons have a stronger H2Bubl signal than flanl<ing exons. (A-F) H2Bubl , H3K36me3, and H3K79me2 signals in skipped exons and 
their flanking exons. (Gray) Three constitutive exons used as controls. In the control, H2Bubl shows decreased patterns from 5' to 3', except in Drosophila 
(n = 16,724 in NCCIT). (Red, black, green) Skipped exon, upstream of the skipped exon, and downstream from the skipped exon, respectively, (/-axis) 
Average modification signals with standard deviation error bars. In NCCIT (n = 1 2,337), HeLa (n = 9722), and myoblast (n = 1 2,356) cells, the H2Bub1 
signals were higher in skipped exons than in flanking exons. In Drosophila, however, the H2Bubl signal was lower in skipped exons than in flanking exons 
(n = 3201 ). The H3K36me3 signal was lower in skipped exons and H3K79me2 did not show any specific pattern in the skipped exons. (C) Differences in 
H2Bub1 between wild-type and siRNF20-transfected NCCIT cells for persistently skipped exons and included exons with standard deviation error bars 
(left). H2Bub1 enrichment is shown for both persistently skipped exons and included exons in siRFN20-transfected cells with standard deviation error bars 
(right). H2Bubl enrichment was determined by dividing the H2Bubl tag density by the average tag density across the whole genome (see Methods). 
P-values were calculated using the KS-test. 
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knockdown of RNF20 expression elicited 
changes to 37% of the skipped exons in 
wild-type cells, and this effect may be at- 
tributable to a reduction of H2Bubl signals. 
When we compared H2Bubl signals be- 
tween wild-type and siRNFZO-transfected 
cells, the persistently skipped exons showed 
lower degrees of H2Bubl reduction (Fig. 
4G, left). A similar pattern was observed 
when we compared H2Bubl signals be- 
tween siEGFP-transfected NCCIT cells 
(control) and siRNFZO-trnasfected NCCIT 
cells (Supplemental Fig. 5D). In the 
siRNF20-transfected NCCIT cells, the in- 
cluded exons showed relatively lower 
H2Bubl signal enrichments compared 
with the persistently skipped exons, sug- 
gesting that increased H2Bubl signals are 
correlated with exon skipping (Fig. 4G, 
right). 

To explain the persistent skipping of 
exons despite reductions in the H2Bubl 
signal, we hypothesized that other mod- 
ifications, such as nucleosome occupancy 
and H3K36me3, may compensate in exon 
skipping. The patterns of H2Bubl and 
nucleosome in wild-type cells seemed 

to be independent of whether the exons were persistently skipped 
or not (Supplemental Fig. 5F,G). However, the H3K36me3 signals 
seemed to correlate more strongly with the persistently skipped 
exons (green lines in Supplemental Fig. 5E, 15% decreased H3K36 
signals compared with flanking exons) compared with included 
exons (blue lines in Supplemental Fig. 5E, 8% decreased H3K36me3 
signals compared with flanking exons), suggesting that H3K36me3 
may compensate for the reduction of H2Bubl. Therefore, various 
chromatin features should be considered together when we seek 
to elucidate the effect of H2Bubl on exon skipping. Although 
there are many limitations to this analysis, our results strongly 
support the idea that H2Bubl plays a potentiating role in exon 
skipping. 

Specific H2Bubl properties in mammals are associated 
with exon skipping 

The correlative patterns between H2Bubl signals and exon skip- 
ping raised many interesting questions. To examine the discrep- 
ancy of H2Bubl patterns in skipped exons between mammals and 
Drosophila, we investigated H2Bubl patterns in ORF regions in hu- 
man, mouse, and fly cells. In NCCIT cells, H2Bubl and H3K36me3 
demonstrated opposing patterns, in that H2Bubl was biased to the 
5' end and H3K36me3 was biased to the 3' end (Fig. 5 A). A similar 
pattern was observed in mouse cells (Fig. 5B). In Drosophila cells, 
however, H3K36me3 demonstrated a pattern similar to that in hu- 
man and mouse cells (Can et al. 2010), while H2Bubl was enriched 
in central ORF regions (Fig. 5C; Kharchenko et al. 2010). The cen- 
tral enrichment of H2Bubl has also been reported in Arabidopsis 
(Roudier et al. 201 1). This phenomenon may suggest that the 5 ' bias 
of H2Bubl in mammals results from a gain-of-function during 
evolution. In light of the observed discrepancy in exon skipping 
between mammals and fly, we suspect that different mechanisms 
underlie the enrichment of H2Bubl in ORF regions, leading to 
different correlation patterns with exon skipping. 
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With regard to the opposing correlative patterns between 
H2Bubl and H3K36me3 in skipped exons, we hypothesized that 
these may arise through differences in the regulation of the nu- 
cleosome architecture. A recent study reported that H2Bubl dis- 
rupts local and higher-order chromatin compaction, indicating 
that it plays significant roles in nucleosome architecture (Fierz 
et al. 2011). In yeast, a functional interplay between SPT16 and 
H2Bubl regulates nucleosome reassembly in GALl gene (Fleming 
et al. 2008), and dynamic regulation of H2B ubiquitylation/ 
deubiquitylation could be equally important for transcription 
(Henry et al. 2003; Daniel et al. 2004; Weake and Workman 2008; 
Song and Ahn 2010). Notably, our data indicate that the H2Bubl 
peak-containing regions demonstrated reduced nucleosome oc- 
cupancy compared with H3K36me3 (Fig. 5D; see Supplemental 
Methods). The mouse genome also demonstrated a similar pattern 
(Fig. 5E). Thus, in the wake of transcription, H2Bubl and H3K36me3 
may affect nucleosome dynamics in different ways, potentially 
leading to opposing correlative patterns of exon skipping. 

The enrichment of H2Bubl in skipped exons (Fig. 4A,B,E) 
seemed to be counterintuitive, given that the H2Bubl signals were 
enriched in highly expressed exons (Fig. 3B,C). A specific set of 
histone modifications (including H2Bubl) may cause a well-de- 
fined "exonic nucleosome" during transcription elongation (Keren 
et al. 2010). The more frequent occurrence of exon skipping in 
highly expressed genes (Fig. 5F) may be one reason for the seem- 
ingly counterintuitive pattern observed herein. We expect that the 
more complex functions of H2Bubl may be linked to exon skip- 
ping in a context-dependent manner. 

Discussion 

The findings in this study suggest that H2B monoubiquitylation plays 
several roles during transcription, acting in both H3 methylati on- 
dependent and -independent fashions. First, H2Bubl is a 5 '-enriched 
active transcription mark in the human genome. Second, H2Bubl 
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has genome-wide correlation patterns with the H3K79 methylations 
and co-occupancies between H2Bubl, and the methylations of 
H3K79 are associated with highly expressed genes. Third, H2Bubl is 
strongly enriched at the exon boundaries of highly expressed exons. 
Fourth, H2Bubl displays an increased density in skipped exons. 
Together, these results define additional functions for H2B mono- 
ubiquitylation during active transcription in humans. 

One interesting finding of this study is that combinatorial 
patterns of H2Bubl and H3K79 methylations are characterized 
by distinct biological functions. Although co-occupancies between 
H2Bubl and the methylations of H3K79 seemed to be strongly 
associated with pre-mRNA processing-related genes (Supplemental 
Fig. lA), the functions of these genes differed according to the 
combinatorial patterns (Fig. 2D). With respect to gene expression 
levels, H3K79me2 showed a stronger tendency for cooperation 
with H2Bubl compared with H3K79 mono- and trimethylations 
(Fig. 2B; Supplemental Fig. 3F,G). Although DOTIL has been 
established to modify all H3K79 methylations, our results suggest 
that the multiple methylations of H3K79 may play distinct roles in 
different contexts. 

In terms of genome-wide occupancy, the ability of H2Bubl to 
stimulate H3K79 methylations is evolutionarily conserved from 
yeast to human. Therefore, like H2Bubl, the methylations of H3K79 
are expected to be enriched surrounding the intron-exon bound- 
aries. However, our present results suggest that H2Bubl has a more 
unique role in determining exon-intron structure. A recent study 
suggested that the effects of H2Bubl in nucleosome stability and 
structure may be related to the enrichment of H2Bubl in intron- 
exon boundaries (Chandrasekharan et al. 2009). Given all of the 
available evidence, the present study sheds light on the role of 
H2Bubl in exon-intron structure, suggesting that it helps medi- 
ate complex splicing mechanisms in higher eukaryotes, inde- 
pendent of the crosstalk between H2Bubl and the methylations 
of H3K79. 

Two models explain the mechanism of exon and intron se- 
lection: intron definition and exon definition (Keren et al. 2010). 
Placing the basal splicing machinery across introns is regarded 
as the ancient mechanism. In mammals, however, although the 
length of the exons has remained relatively constant, intron length 
has gradually expanded. As opposed to fungi, in which introns have 
remained short and are thought to be the recognized unit for se- 
lection (i.e., intron definition), higher eukaryotes use exon defi- 
nition as their major mechanism for selection. In this regard, H2Bubl 
can be hypothesized to have achieved an uncharacterized function 
in exon definition due to evolutionary pressure. Further studies 
will be required to validate our hypothesis through comparative 
genomic approaches. 

Ubiquitylation, which results in the addition of a bulky 
moiety to lysine residues within proteins, exerts a dramatic effect 
on protein conformation and function. Monoubiquitylation of 
the H2B C terminus affects chromatin structure in vivo, thereby 
modulating transcription elongation (Weake and Workman 2008). 
Here, we found that histone H2Bubl is enriched at exon boundaries 
(Fig. 3; Supplemental Fig. 4) and correlates with exon skipping 
(Fig. 4). According to the "speed bump" model, the speed of RNA 
PolII determines the binding of splicing factors during transcrip- 
tion. In this model, the speed of PolII is regulated by various barriers, 
including high-nucleosome density, DNA methylation, and several 
types of histone modifications (Keren et al. 2010; Schwartz and Ast 
2010). All of the present correlative findings reinforce the notion 
that H2Bubl potentiates alternative splicing. Moreover, we found 
that the methylations of H3K4 and H3K79 are not tightly related to 



exon-intron structure or exon skipping in human cells, suggesting 
that H2Bubl plays a unique role in chromatin architecture, inde- 
pendent of the methylations of histone H3. 

Methods 

Antibodies 

Polyclonal rabbit antibodies recognizing H3ac, H3K4me3, 
H3K36me3, H3K79mel, H3K79me2, and H3K79me3 were pur- 
chased from Abeam or produced in-house as previously described 
(Oh et al. 2010). Antibodies against H2Bubl were obtained from 
Millipore (17-650) and Cell Signaling (5546). The antibody speci- 
ficity of the anti-H2Bubl from Millipore is shown in Supplemental 
Figure 6. The anti-H2Bubl from Cell Signaling was also confirmed 
(data not shown). 

Cell culture, chromatin immunoprecipitation [CHIP) assay, 
and mononucleosome purification 

The NCCIT embryonic carcinoma cell line was cultured on 10-cm 
dishes in RPMI 1640 medium supplemented with 10% fetal bovine 
serum (FBS, Hy clone) and 1% antibiotic-antimycotic solution 
(Hyclone) at 37°C in a humidified atmosphere containing of 5% CO2. 

Cross-linking was performed by incubating cells in 1% form- 
aldehyde (Sigma) for 20 min at 25°C and then quenching cells with 
125 mM glycine. For ChIP, phosphate buffered saline (PBS) -washed 
cells were lysed in 400 jxL of SDS lysis buffer (1% SDS, 10 mM EDTA, 
50 mM Tris-HCl at pH 8.0), and DNA was sonicated to create frag- 
ments 400-500 nt in length using a Sonic Dismembrator (Model 
500, Fisher). The chromatin was cleared by centrifugation for 15 min 
at 13,000 rpm and 4°C and subjected to immunoprecipitation for 
2 h at 4°C using the indicated antibodies and a mixture of protein 
A and G agarose (GE Healthcare). The protein A and G agarose/ 
antibody/ chromatin complexes were washed for 10 min each with 
the following solutions: low-salt wash buffer (20 mM Tris-Cl at pH 
8.0, 0.1% SDS, 2 mM EDTA, 1% Triton X- 100, ISOmMNaCl), high- 
salt wash buffer (20 mM Tris-Cl at pH 8.0, 0.1% SDS, 2 mM EDTA, 
1% Triton X-100, 500 mM NaCl), LiCl buffer (10 mM Tris-Cl at pH 
8.0, 250 mM LiCl, 1 mM EDTA, 1% NP-40, 1% deoxycholic acid), 
and TE buffer (10 mM Tris-Cl at pH 8.0, 1 mM EDTA). The chro- 
matin complexes were eluted twice in 250 jxL of elution buffer (1% 
SDS, 0.1 M NaHCOa) by incubation at 25°C for 15 min with ro- 
tation, and then reverse cross-linked by heating overnight at 68°C 
with 400 mM NaCl. 

For mononucleosome purification, PBS-washed cells were 
lysed in 400 fjiL of lysis buffer (18% Ficoll-400 [Sigma F2637], 
10 mM KH2PO4, 10 mM K2PO4, 1 mM MgCl2, 0.25 mM EGTA, and 
0.1 mM PMSF) and centrifuged at 4°C for 15 min at 13,000 rpm. 
The nuclear pellets were resuspended with 1 mL of Buffer A (10 mM 
Tris-Cl at pH 7.4, 150 mM NaCl, 5 mM KCl, 1 mM EDTA, 0.1 mM 
PMSF), and then mixed with 5 |jlL of CaCl2 (final concentration, 
5 mM) and 1.5 |jlL of MNase (500 U/2.5 |jlL, NEB). After incubation 
at 37°C for 5 min, the reactions were stopped by adding 100 |jlL of 
1% SDS and 50 |jlL of 0.5 mM EDTA. The MNase-digested nuclei 
were reverse cross-linked by heating overnight at 68°C with 350 mM 
NaCl. Finally, the reverse cross-linked samples were treated with 
RNase A and proteinase K (Sigma), isolated by phenol/chloroform 
extraction, and precipitated using 3 M sodium acetate (pH 5.2) with 
absolute ethanol. 

slRNA and transfection 

A 21-nt synthetic siRNA duplex designed to inhibit expression of 
RNF20 (sense, 5'-ACAUCCGUAUCAUCCUUAAUU-3', and anti- 
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sense, 5'-UUAAGGAUGAUACGGAUGUUU-3') was synthesized 
and purified by Samchully Pharm. The day before transfection, 
NCCIT cells were seeded at a density of 2 x 10^ cells/well in 6-well 
plates. siRNFZO was then transfected into the cells using FuGene 
HD (Roche) according to the manufacturer's instructions. Three days 
after transfection, the samples were separated by SDS-PAGE and 
subjected to immunoblot analysis using the indicated antibodies. 

Illumina library construction and sequencing 

The purified ChIP DNA fragments were ligated to a pair of adaptors 
designed to enable sequencing on the Illumina Genome Analyzer. 
The ligation products were size-fractionated on an agarose gel, and 
the 200-300-bp fragments were subjected to 18 cycles of PGR 
amplification. Cluster generation and 72 cycles of single-read 
sequencing were performed. 

The mRNA sequencing library was generated using the mRNA- 
seq sample prep kit (Illumina). Polyadenylated mRNA was purified 
from 10 fxg of total RNA using Sera-mag oligo(dT) beads (Illumina), 
and then fragmented by incubation with fragmentation buffer 
(Illumina) at 94°C for 5 min. First-strand cDNA was synthesized 
using Superscript II (Invitrogen) and random primers (Illumina), 
and second-strand cDNA was synthesized using RNaseH and DNA 
polymerase I (Illumina). The resulting cDNA fragments were ligated 
to a pair of adaptors designed for sequencing on the Illumina Ge- 
nome Analyzer. The ligation products were size-fractionated to ob- 
tain 200-bp fragments on a 2% agarose gel and subjected to 15 cycles 
of PGR amplification. Cluster generation and 72 cycles of single-read 
sequencing were performed. 

Data processing 

The sequence tags were mapped to the human genome (UCSC 
hgl8) using Bowtie (Langmead et al. 2009) with default options. 
The mapped sequence reads were extended to the average fragment 
length (200 bp for ChlP-seq and 150 bp for MNase-seq), and the 
number of overlapping sequence tags was determined using 30-bp 
windows with 15-bp steps across the human genome. The ChlP-seq 
signal was defined as the ratio of target read count/target size di- 
vided by total read count/genome size Qung and Kim 2009). The 
mRNA-seq data were mapped to the human genome using TopHat 
(Trapnell et al. 2009) with default parameters, and the expression 
level of each gene was defined using the FPKM (fragment per kilo- 
base of exon per million fragments mapped) scores generated by 
Cufflinks (Trapnell et al. 2010). Genes with FPKM scores <0.05 were 
removed, as they could represent unexpressed or unsequenced re- 
gions. The expression level of each exon was defined as exon read 
count/exon size divided by total read count/genome size. 

Data analysis 

Genome annotations were downloaded from the UCSC Genome 
Browser, human genome Build 36.1 (hglS assembly). Gene defi- 
nitions were obtained from Refseq Gene track in the UCSC Ge- 
nome Browser (Karolchik et al. 2009). Exons were defined using 
Refseq Gene track and only internal exons separated from flanking 
exons by at least 500 bp were analyzed. 

For HeLa cells, the H2Bubl ChlP-seq data were kindly provided 
by Shema et al. (2008). The mRNA-seq data were downloaded from the 
GEO database with accession number GSM546927 (Guo et al. 2010). 

For our analysis of Drosophila melanogaster, H2Bubl ChlP-chip 
and mRNA-seq for the S2-DRSC cell line (late embryonic stage em- 
bryo-derived cell line) were downloaded from the GEO database 
using accession numbers GSE20773 (Roy et al. 2011) and 
GSM461176 (Brooks et al. 2011), respectively. The H3K36me3 



ChlP-seq data were downloaded from the GEO database with ac- 
cession number GSM480158 (Gan et al. 2010). Gene definitions 
were given by their Refseq gene track (dm3) in the UCSC Genome 
Browser and other processes were the same as those used for the 
human genome data. 

For mouse myoblast cells, the H2Bubl and H3K36me3 ChlP- 
seq data were downloaded from ArrayExpress with accession num- 
ber E-GEOD-25308 (Asp et al. 2011; Vethantham et al. 2012), and 
the RNA-seq data were downloaded from the GEO database with 
accession number GSM521256 (Trapnell et al. 2010). Gene defini- 
tions were given by their Refseq gene track (mm9) in the UCSC 
Genome Browser, and the other processes were the same as those 
used for the human genome data. 

Calculating modification signals in skipped exons 

Alternatively spliced exons were defined by more than threefold 
repression of cassette exons compared with flanking exons after 
the removal of the first and last two exons from the analysis. To 
precisely define skipped exons, only exons of highly expressed 
genes (the top 30%) were considered; however, the same tendency 
was achieved when all genes were used. In NCCIT cells, HeLa 
cells, siRNF20-transfected NCCIT cells, mouse myoblast cells, and 
Drosophila melanogaster cells, we observed 12337, 9722, 11496, 
12356, and 3201 skipped exons, respectively. As a control, we se- 
lected three constitutive exons that were separated from skipped 
exons by at least two exons. NCCIT cells contained 16,724 control 
exons. All exon definitions were provided by the Refseq gene tracks 
for all organisms, as obtained from the UCSC Genome Browser. 

When we compared exon skipping between wild-type and 
siRNF20-transfected cells, we selected skipped exons with H2Bubl- 
specific patterns in order to infer the effect of RNF20 knockdown. 
Among the skipped exons, we selected those showing a more than 
or less than 1.5-fold enrichment of H2Bubl signals compared with 
flanking exons. We found that 5821 exons were persistently skipped 
in both cell types, while 2749 exons were skipped only in wild-type 
cells. To check the effect of our definition of exon skipping, we 
changed the threshold value for exon expression levels from 
threefold to twofold. Using this parameter, we observed a similar 
tendency in our results. The H2Bubl signal difference between 
wild-type and siRNF20-transfected cells was calculated by com- 
paring the fold enrichment of H2Bubl signals versus the average 
H2Bubl signal (total read count/genome size). 

Data access 

The ChlP-seq data from this study have been submitted to the NCBI 
Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/ 
geo/) under accession number GSE25882. 

Acknowledgments 

Inkyung Jung thanks the members of the Bioinformatics and 
Computational Biology Laboratory and Yoo-Sun Lee for helpful 
discussions. We also thank the modENCODE project for allowing 
us to use their data. This work is also supported by grants from 
the Epigenomic Research Program for Human Stem Cells (2007- 
2004134), the Future-Based Technology Development Program of 
the National Research Foundation of Korea (2010-0020590), and 
the KAIST Future Systems Healthcare Project, Ministry of Educa- 
tion, Science and Technology, South Korea. This work is also sup- 
ported by the CHUNG Moon Soul Center for Biolnformation and 
BioElectronics (CMSC), and by the Ministry of Education, Science 
and Technology, South Korea (Bioinformatics Pipeline for Stem 
Cell Epigenomics). 



1 034 Genome Research 

www.genome.org 



H2Bubl determines intron-exon boundaries 



References 

Andersson R, Enroth S, Rada-Iglesias A, Wadelius C, Komorowski J. 2009. 

Nucleosomes are well positioned in exons and carry characteristic 

histone modifications. Genome Res 19: 1732-1741. 
Asp P, Blum R, Vethantham V, Parisi F, Micsinai M, Cheng J, Bowman C, 

Kluger Y, Dynlacht BD. 2011. PNAS Plus: Genome-wide remodeling of 

the epigenetic landscape during myogenic differentiation. Proc Natl 

Acad Sci 108: E149-E158. 
Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev 

1; Zhao K. 2007. High-resolution profiling of histone methylations in the 

human genome. Cell 129: 823-837. 
Batta K, Zhang Z, Yen K, Goffman DB, Pugh BF. 2011. Genome-wide 

function of H2B ubiquitylation in promoter and genie regions. Genes 

Dev25: 2254-2265. 
Briggs SD, Xiao X Sun ZW, Caldwell JA, Shabanowitz J, Hunt DF, AUis CD, 

Strahl BD. 2002. Gene silencing: trans-histone regulatory pathway in 

chromatin. Nature ^IS: 498. doi: 10.1038/nature00970. 
Brooks AN, Yang L, Duff MO, Hansen KD, Park JW, Dudoit S, Brenner SE, 

Graveley BR. 2011. Conservation of an RNA regulatory map between 

Drosophila and mammals. Genome Res 21: 193-202. 
Chandrasekharan MB, Huang F, Sun ZW. 2009. Ubiquitination of histone 

H2B regulates chromatin dynamics by enhancing nucleosome stability. 

Proc Natl Acad Sci 106: 16686-16691. 
Choi JK. 2010. Contrasting chromatin organization of CpG islands and 

exons in the human genome. Genome Biol 11: R70. doi: 10.1186/gb- 

2010-ll-7-r70. 

Damjanov 1. 1993. Pathogenesis of testicular germ cell tumours. Eur Urol 23: 
2-5. 

Damjanov 1, Horvat B, Gibas Z. 1993. Retinoic acid-induced differentiation 
of the developmentally pluripotent human germ cell tumor-derived cell 
line, NCCIT. Lab Invest 68: 220-232. 

Daniel JA, Torok MS, Sun ZW, Schieltz D, Allis CD, Yates JR 3rd, Grant PA. 
2004. Deubiquitination of histone H2B by a yeast acetyltransferase 
complex regulates transcription. J Biol Chem 279: 1867-1871. 

Dover J, Schneider J, Tawiah-Boateng MA, Wood A, Dean K, Johnston M, 
Shilatifard A. 2002. Methylation of histone H3 by COMPASS requires 
ubiquitination of histone H2B by Rad6. J Biol Chem 277: 28368-28371. 

Fierz B, Chatter] ee C, McGinty RK, Bar-Dagan M, Raleigh DP, Muir TW. 
2011. Histone H2B ubiquitylation disrupts local and higher-order 
chromatin compaction. Nat Chem Biol 7: 113-119. 

Fleming AB, Kao CF, Hillyer C, Pikaart M, Osley MA. 2008. H2B 
ubiquitylation plays a role in nucleosome dynamics during 
transcription elongation. Mol Cell 31: 57-66. 

Gan Q, Schones DE, Ho Eun S, Wei G, Cui K, Zhao K, Chen X. 2010. 

Monovalent and unpoised status of most genes in undifferentiated cell- 
enriched Drosop/zi/a testis. Genome Biol 11: R42. doi: 10.1186/gb-2010- 
ll-4-r42. 

Guccione E, Bassi C, Casadio F, Martinato F, Cesaroni M, Schuchlautz H, 
Luscher B, Amati B. 2007. Methylation of histone H3R2 by PRMT6 and 
H3K4 by an MLL complex are mutually exclusive. Nature 449: 933-937. 

Guo H, Ingolia NT, Weissman JS, Bartel DP. 2010. Mammalian micro RNAs 
predominantly act to decrease target mRNA levels. Nature 466: 835- 
840. 

Henry KW, Wyce A, Lo WS, Duggan LJ, Emre NC, Kao CF, Pillus L, Shilatifard 

A, Osley MA, Berger SL. 2003. Transcriptional activation via sequential 

histone H2B ubiquitylation and deubiquitylation, mediated by SAGA- 

associated Ubp8. Genes Dev 17: 2648-2663. 
Hon G, Wang W, Ren B. 2009. Discovery and annotation of functional 

chromatin signatures in the human genome. PLoS Comput Biol 5: 

el000566. doi: 10.1371/journal.pcbi.l000566. 
Hung T, Binda O, Champagne KS, Kuo AJ, Johnson K, Chang HY, Simon 

MD, Kutateladze TG, Gozani O. 2009. 1NG4 mediates crosstalk between 

histone H3 K4 trimethylation and H3 acetylation to attenuate cellular 

transformation. Mol Cell 33: 248-256. 
Jung 1, Kim D. 2009. Regulatory patterns of histone modifications to control 

the DNA methylation status at CpG islands. IBC 1: 1-8. 
Karolchik D, Hinrichs AS, Kent WJ. 2009. The UCSC Genome Browser. Curr 

Protoc Bioinformatics 28: 1.4.1-1.4.26. 
Keren H, Lev-Maor G, Ast G. 2010. Alternative splicing and evolution: 

Diversification, exon definition and function. Nat Rev Genet 11: 345-355. 
Kharchenko PV, Alekseyenko AA, Schwartz YB, Minoda A, Riddle NC, Ernst 

J, Sabo PJ, Larschan E, Gorchakov AA, Gu T, et al. 2010. Comprehensive 

analysis of the chromatin landscape in Drosophila melanogaster. Nature 

471: 480-485. 

Kim J, Hake SB, Roeder RG. 2005. The human homolog of yeast BREl 
functions as a transcriptional coactivator through direct activator 
interactions. Mol Cell 20: 759-770. 

Kim J, Guermah M, McGinty RK, Lee JS, Tang Z, Milne TA, Shilatifard A, 
Muir TW, Roeder RG. 2009. RAD6-Mediated transcription-coupled H2B 



ubiquitylation directly stimulates H3K4 methylation in human cells. 
Cell 137: 459-471. 

Kolasinska-Zwierz P, Down T, Latorre 1, Liu T, Liu XS, Ahringer J. 2009. 
Differential chromatin marking of introns and expressed exons by 
H3K36me3. Nat Genet 41: 376-381. 

Langmead B, Trapnell C, Pop M, Salzberg SL. 2009. Ultrafast and memory- 
efficient alignment of short DNA sequences to the human genome. 
Genome Biol 10: R25. doi: 10.1186/gb-2009-10-3-r25. 

Luco RF, Pan Q, Tominaga K, Blencowe BJ, Pereira-Smith OM, Misteli T. 
2010. Regulation of alternative splicing by histone modifications. 
Science 327: 996-1000. 

McGinty RK, Kim J, Chatterjee C, Roeder RG, Muir TW. 2008. Chemically 
ubiquitylated histone H2B stimulates hDotlL-mediated 
intranucleosomal methylation. Nature 453: 812-816. 

Minsky N, Shema E, Field Y, Schuster M, Segal E, Oren M. 2008. 

Monoubiquitinated H2B is associated with the transcribed region of 
highly expressed genes in human cells. Nat Cell Biol 10: 483-488. 

Oh S, Jeong K, Kim H, Kwon CS, Lee D. 2010. A lysine-rich region in Dotlp 
is crucial for direct interaction with H2B ubiquitylation and high 
level methylation of H3K79. Biochem Biophys Res Commun 399: 512- 
517. 

Pavri R, Zhu B, Li G, Trojer P, Mandal S, Shilatifard A, Reinberg D. 2006. 
Histone H2B monoubiquitination functions cooperatively with FACT to 
regulate elongation by RNA polymerase 11. Cell 125: 703-717. 

Roudier F, Ahmed 1, Berard C, Sarazin A, Mary-Huard T, Cortijo S, Bouyer D, 
Caillieux E, Duvernois-Berthet E, Al-Shikhley L, et al. 2011. Integrative 
epigenomic mapping defines four main chromatin states in Arabidopsis. 
EMBOf 30: 1928-1938. 

Roy S, Ernst J, Kharchenko PV, Kheradpour P, Negre N, Eaton ML, Landolin 
JM, Bristow CA, Ma L, Lin MF, et al. 2011. Identification of functional 
elements and regulatory circuits by Drosophila modENCODE. Science 
330: 1787-1797. 

Schulze JM, Hentrich T, Nakanishi S, Gupta A, Emberly E, Shilatifard A, 
Kobor MS. 2011. Splitting the task: Ubp8 and UbplO deubiquitinate 
different cellular pools of H2BK123. Genes Dev 25: 2242-2247. 

Schwartz S, Ast G. 2010. Chromatin density and splicing destiny: On the 
cross-talk between chromatin structure and splicing. EMBO J 29: 1629- 
1636. 

Shahbazian MD, Zhang K, Grunstein M. 2005. Histone H2B ubiquitylation 
controls processive methylation but not monomethylation by Dotl and 
Setl. Mo/ Cell 19: 271-277. 

Shema E, Tirosh 1, Aylon Y, Huang J, Ye C, Moskovits N, Raver-Shapira N, 
Minsky N, Pirngruber J, Tarcic G, et al. 2008. The histone H2B-specific 
ubiquitin Ugase RNF20/hBREl acts as a putative tumor suppressor through 
selective regulation of gene expression. Genes Dev 22: 2664-2676. 

Song YH, Ahn SH. 2010. A Brel -associated protein, large 1 (Lgel), promotes 
H2B ubiquitylation during the early stages of transcription elongation. 
J Biol Chem 285: 2361-2367. 

Stephens MA. 1970. Use of the Kolmogorov-Smirnov, Cramer- Von Mises 
and related statistics without extensive tables. JR StatSoc 32: 115-122. 

Sun ZW, Allis CD. 2002. Ubiquitination of histone H2B regulates H3 
methylation and gene silencing in yeast. Nature 418: 104-108. 

Tanny JC, Erdjument-Bromage H, Tempst P, Allis CD. 2007. Ubiquitylation 
of histone H2B controls RNA polymerase 11 transcription elongation 
independently of histone H3 methylation. Genes Dev 21: 835-847. 

Teshima T, Chatani M, Hata K, Inoue T. 1988. High-dose rate intracavitary 
therapy for carcinoma of the uterine cervix: 11. Risk factors for rectal 
complication. IntJRadiat Oncol Biol Phys 14: 281-286. 

Tilgner H, Nikolaou C, Althammer S, Sammeth M, Beato M, Valcarcel J, 
Guigo R. 2009. Nucleosome positioning as a determinant of exon 
recognition. Nat Struct Mol Biol 16: 996-1001. 

Trapnell C, Pachter L, Salzberg SL. 2009. TopHat: Discovering splice 
junctions with RNA-Seq. Bioinformatics 25: 1105-1111. 

Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, 
Salzberg SL, Wold BJ, Pachter L. 2010. Transcript assembly and 
quantification by RNA-Seq reveals unannotated transcripts and isoform 
switching during cell differentiation. Nat Biotechnol 28: 511-515. 

Vethantham V, Yang Y, Bowman C, Asp P, Lee JH, Skalnik DC, Dynlacht BD. 
2012. Dynamic loss of H2B ubiquitylation without corresponding 
changes in H3K4 trimethylation during myogenic differentiation. Mol 
Cell Biol 32: 1044-1055. 

Wang Z, Zang C, Rosenfeld JA, Schones DE, Barski A, Cuddapah S, Cui K, 
Roh TY, Peng W, Zhang MQ, et al. 2008. Combinatorial patterns of 
histone acetylations and methylations in the human genome. Nat Genet 
40: 897-903. 

Weake VM, Workman JL. 2008. Histone ubiquitination: Triggering gene 
activity. Mol Cell 29: 653-663. 



Received January 9, 201 1; accepted in revised form March 7, 2012. 



Genome Research 1 035 

www.genonne.org 



